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The dark matter interpretation of the DAM A modulation signal depends on the Nal(Tl) scintil¬ 
lation efficiency of nuclear recoils. Previous measurements for Na recoils have large discrepancies, 
especially in the DAMA/LIBRA modulation energy region. We report a quenching effect measure¬ 
ment of Na recoils in Nal(Tl) from 3keV nr to 52keV nr , covering the whole DAMA/LIBRA energy 
region for light WIMP interpretations. By using a low-energy, pulsed neutron beam, a double 
time-of-flight technique, and pulse-shape discrimination methods, we obtained the most accurate 
measurement of this kind for Nal(Tl) to date. The results differ significantly from the DAMA 
reported values at low energies, but fall between the other previous measurements. We present 
the implications of the new quenching results for the dark matter interpretation of the DAMA 
modulation signal. 


I. INTRODUCTION 

For over a decade, the DAMA experiments (DAMA- 
Nal and DAMA/LIBRA) have been observing an annual 
modulation in the rate of events in the low-energy region 
of Nal(Tl) detectors pQ. This modulation signal has an 
extremely high statistical significance (9.3cr) and is often 
interpreted as evidence for WIMP dark matter interac¬ 
tions, such as low mass WIMP scattering [2] or inelastic 
WIMP scattering 011 ]. Several experiments have ruled 
out the DAMA dark matter claim in the standard WIMP 
picture [SHE], while alternative WIMP theories might still 
reconcile the experimental results mm- 

The dark matter interpretation of the DAMA mod¬ 
ulation signal depends on the scintillation efficiency of 
Nal(Tl) for sodium and iodine recoils relative to that 
of gammas (electron recoils); the former could be in¬ 
duced by WIMP scattering interactions, while the lat¬ 
ter are used to calibrate the detectors. Nuclear recoils 
in Nal(Tl), due to the small fraction of energy trans¬ 
fer to electrons, typically produce less scintillation light 
compared with electron recoils with the same energy de¬ 
position. The relative ratio is usually referred to as the 
nuclear recoil quenching factor. This factor is critical 
to translate the observed electron equivalent energy into 
nuclear recoil energy for dark matter analysis. 

DAMA reports a quenching factor of 0.3 for sodium 
recoils and 0.09 for iodine recoils m • These results were 
obtained by exposing a Nal(Tl) detector to neutrons from 
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a 252 Cf source. The nuclear recoil signals produced by 
the neutrons were compared to Monte Carlo-simulated 
sodium and iodine recoil spectra to extract the quenching 
factors, which were assumed to be energy-independent. 
Since then, several more experiments have been carried 
out to measure the sodium and iodine quenching fac¬ 
tors as a function of nuclear recoil energy using mono- 
energetic neutron sources, mostly deuterium-deuterium 
neutron generators [121416] . By looking for coincidence 
signals between a Nal(Tl) detector and neutron detec¬ 
tors at fixed neutron scattering angles, this technique 
could provide a direct measurement of the DAMA mod¬ 
ulation energy scale if the signal is interpreted as nuclear 
recoils. A few of these measurements reported quench¬ 
ing factors consistent with the DAMA results, especially 
in energy regions higher than 20keV nuclear recoil en¬ 
ergy (20keV nr ). However, recent measurements by Cha- 
gani m and Collar [18] led to new Na recoil quench¬ 
ing factors significantly deviating from the DAMA val¬ 
ues over a wide energy range, and these two new mea¬ 
surements also conflict seriously with each other at low 
energies, where the DAMA modulation signals occur in 
a light WIMP interpretation. 

In addition to these inconsistencies, the previous Na 
recoil quenching measurements in Nal(Tl) typically carry 
large uncertainties, ranging from ~ 10% (relative) to, at 
the lowest energies, over 100% (relative). Therefore, it is 
necessary to conduct new quenching measurements that 
can significantly improve the quenching-factor accuracy 
and provide a reliable Na recoil calibration for the light 
WIMP interpretation of the DAMA results, as well as for 
other Nal(Tl) dark matter experiments [191423] . 

In this paper, we report a Nal(Tl) quenching measure- 
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ment using a pulsed neutron beam produced by the FN 
tandem facility at the University of Notre Dame Nuclear 
Science Laboratory. This measurement was designed to 
achieve an overall uncertainty of ^5% and an energy 
threshold of a few keV nr . Several techniques were com¬ 
bined to suppress backgrounds and uncertainties, as sum¬ 
marized below: 

1. A triple time-coincidence between a pulsed neutron 
beam, a Nal(Tl) detector and an angular array of 
neutron detectors was used to select neutron events 
and to reduce random coincidence backgrounds. 

2. Low-energy neutrons 690 keV) were used so that 
low-energy nuclear recoils (< 50 keV nr ) could be ob¬ 
tained at large neutron scattering angles and the 
relative angular uncertainties were reduced. 

3. A small Nal(Tl) crystal (25 mm cube) was used to 
reduce multiple scattering backgrounds. 
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4. A high-quant urn-efficiency photomultiplier 

(PMT) was used to enhance light collection 

18photoelectrons/keV ee achieved). 

5. Pulse-shape discrimination (PSD) methods were 
used to select neutrons events and to reject gammas 
and noise. 

This experiment was inspired by the success of the 
SCENE experiments, which measured the nuclear recoil 
quenching effects in liquid argon down to very low nu¬ 
clear recoil energy (^10keV nr ) using the Notre Dame 
facility |24j [25] . 


FIG. 1. To-scale, bird’s-eye view of the experimental setup in 
the first position configuration with side view of the Nal(Tl) 
detector. A LiF target in the beam produces neutrons when 
struck by the proton beam. The neutrons travel to the 
Nal(Tl) detector, where they scatter off a sodium or iodine 
nucleus to one of the neutron detectors positioned at the rel¬ 
evant angle. Two different types of neutron detectors were 
used, as described in the text below. A polyethylene colli¬ 
mator prevents neutrons from hitting the neutron detectors 
directly from the LiF target. Also shown is a side view of 
the Nal(Tl) detector, with the crystal shown in white and the 
PMT in green. 


II. EXPERIMENTAL SETUP 
A. Overview 


measuring 12 nuclear recoil energies. A scheme of the 
experimental setup in the first position configuration is 
shown in Figure []] 


The measurement of the Na quenching factor was per¬ 
formed at the FN Tandem accelerator at the University 
of Notre Dame Institute for Structure and Nuclear Astro¬ 
physics. A pulsed beam of protons from the accelerator 
interacted with a LiF target to produce neutrons with a 
nominal energy of 690 keV at 0° scattering angle. A de¬ 
tector consisting of an enclosed Nal(Tl) crystal and pho¬ 
tomultiplier tube (PMT) was placed on the beam line. 
The neutrons traveling in this direction could produce 
nuclear recoil scintillation events in the crystal and be 
subsequently detected by a liquid-scintillator-based neu¬ 
tron detector at a fixed recoil angle. The kinematics of 
this interaction determined the nuclear recoil energy. 

This information, combined with the scintillation light 
collected by the Nal(Tl) detector, provided a measure of 
the light yield of the detector for nuclear recoil events. 
Calibrating this system with electron recoils of known en¬ 
ergies allowed for a determination of the quenching fac¬ 
tor. We used a single Nal(Tl) detector in two positions 
and a stationary array of six neutron detectors, thereby 


B. The Proton Beam and LiF Target 

A beam of 2.44-MeV protons was produced by an 11- 
MeV FN Tandem accelerator. These protons, incident on 
a LiF target, produce neutrons through the 7 Li(p,n) 7 Be 
reaction (Q-value: —1.644 MeV). Ignoring energy loss in 
the target, the neutrons emitted at a nominal scattering 
angle of 0° have an energy of 723 keV. 

The beam was separated into time-bunched pulses by 
a three-part pulsing system with a timing resolution of 
2 ns and an intrinsic period of 101.5 ns. The proton pulse 
selector can additionally be set to allow only one out 
of every n pulses to pass through, effectively increasing 
the period. The beam cross section is 3 mm in diameter 
and is highly stable, with a variation in proton energy of 
around 1 keV. 

The energy of the beam was chosen to increase the 
event rate while reducing the relative angular uncertainty 
for nuclear recoils. Because of the finite detector size of 
both the Nal(Tl) detector and the neutron detectors, as 
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well as the shape of the recoil-energy dependence on an¬ 
gle, the spread in the nuclear recoil energies deposited 
in the crystal is smaller when the neutron detectors are 
located at larger angles with respect to the beam. We 
therefore wanted the proton beam energy to produce the 
recoil energies of interest at relatively large scattering 
angles, while providing enough event rate to collect ade¬ 
quate statistics for the measurement. A calculation was 
done to obtain the overall interaction rate in the Nal(Tl) 
detector at low recoil energies given the 7 Li(p,n) 7 Becross 
sections in [261 and the differential neutron elastic scat¬ 
tering cross sections in 23 Na and 127 I at the appropriate 
angles. A broad maximum in the event rate for all recoil 
energies was found at a proton energy of 2.44 MeV with 
reasonably large scattering angles for recoil energies be¬ 
tween 3 and 52 keV, and therefore that energy was chosen 
for the measurement. 

A higher pulse frequency can lead to a higher neu¬ 
tron event rate, but it may also cause pileups in the Nal 
scintillation time window. Based on a calculation of the 
neutron yield, we determined that a pulse separation of 
609 ns was enough to reduce the pileup rate. As described 
in Section [TTCl the detectors were arrayed in two geomet¬ 
rical configurations. In the first position configuration, 
where the Nal(Tl) detector was 50 cm from the LiF tar¬ 
get, the bunching ratio for the pulser was set to 1 in 6, for 
an effective pulse period of 609 ns, but was later changed 
to 1 in 8 after a high event rate was observed, for a pulse 
period of 812 ns. The second position configuration, with 
a Nal(Tl) distance of 91 cm, had a bunching ratio of 1 in 
8. Each pulse carried ~ 10 4 protons. 

The LiF target was deposited on a 0.4-mm tantalum 
backing, which stops the proton beam and minimizes the 
gamma background. Incoming protons lose energy as 
they travel through the LiF target, leading to a broaden¬ 
ing of the outgoing neutron energy spectrum. The LiF 
target thickness was chosen to be 0.52 mg/cm 2 in order 
to compromise between the event rate and the spread 
in the neutron energy, which both increase with thick¬ 
ness. The mean neutron energy for a target thickness of 
0.52 mg/cm 2 was calculated to be 690 keV with a spread 
of 4% at the full-width-half-maximum. At that thick¬ 
ness, and with a bunching ratio of 1 in 8, the neutron 
flux was calculated to be around 300 neutrons/s at the 
Nal(Tl) detector for the first position configuration and 
about 100 neutrons/s in the second. 


C. The Detectors 

The Nal(Tl) detector consisted of a 25-mm cubical 
Nal(Tl) crystal optically coupled to a 76-mm super- 
bialkali Hamamatsu R6233-100 PMT. The crystal was 
grown at Radiation Monitoring Devices Inc. with high- 
purity Astro-grade Nal powder from Sigma Aldrich. The 
small crystal size was chosen to minimize the probability 
of neutron multiple scattering. The crystal and the PMT 
were packaged in a stainless-steel enclosure with a thin 


wall (^0.5 mm) in the section surrounding the crystal to 
further minimize the chance of multiple scatters. 

A high light-collection efficiency of the detector is nec¬ 
essary to obtain a high energy resolution and low thresh¬ 
old. The PMT had a high peak quantum efficiency of 
^35%. In addition, the crystal was covered on the other 
five faces with highly-reflective Lumirror reflector (>98% 
reflective above 350 nm) additionally wrapped in several 
layers of PTFE tape. No light guide was used in this 
experiment; the coupling was a transparent optical gel 
from Cargille Labs with a refractive index of 1.52. This 
arrangement allowed for a high maximum light yield of 
18.2T0.1 photoelectrons (p.e.)/keV ee , where a keV ee is a 
unit describing the electron-equivalent energy that would 
produce the same scintillation yield as a nuclear recoil 
with an energy of 1 keV nr . 

The Nal(Tl) detector was placed on the beam-line axis 
to maximize the event rate. The Nal(Tl) detector was 
placed 50 cm from the LiF target in the first position 
configuration and 91 cm in the second configuration, as 
described in Table [I| The 50 cm position was chosen in 
order to produce a high event rate while keeping the total 
angular spread, A 6/0, below 5%. The 91-cm position was 
chosen to increase the number of recoil energies explored 
without changing the already well-established locations 
of the neutron detectors. 

The angles chosen for the neutron detectors, between 
18 and 84 degrees, allowed for data to be collected for 
Na nuclear recoil energies between 3 and 52 keV. The 
detector distances were chosen to maximize the event 
rate while maintaining a recoil energy uncertainty due 
to finite detector size of less than 5%. Their positions 
are summarized in Table [TJ Two types of neutron de¬ 
tectors were used for the measurement: 5.1-cm x 05.1- 
cm Eljen 510-20x20-9/301 and 12.7-cm x 012.7-cm El- 
jen 510-50x50-1/301 liquid scintillator detectors. Both 
types have the reflector EJ-510 and the liquid scintillator 
EJ-301, a xylene-based scintillator with organic fluors. 
This scintillator has pulse-shape discrimination capabil¬ 
ity, which allows for the selection of events induced by 
desired particle types. A typical light yield response of 
these detectors was measured to be ^1 p.e./keV ee . 

A 22-cm-diameter, 22-cm-long, cylindrical polyethy¬ 
lene collimator with a 2.5-cm-diameter hole was used to 
prevent neutrons from traveling directly from the LiF 
target to the neutron detectors. 


D. Electronics and Data Acquisition 

The data acquisition system needed to provide an ac¬ 
curate determination of the event energy and timing, as 
well as the particle type through pulse-shape discrimina¬ 
tion. This could be achieved by recording the waveforms 
from the photomultiplier tube signals in both the Nal(Tl) 
and the neutron detectors, as well as the signal from the 
proton pulse selector, during neutron-induced scintilla¬ 
tion events. The data acquisition scheme is shown in 
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TABLE I. Detector information and positions for position 
configurations 1 (top) and 2 (bottom). The “Flight Distance” 
for the neutron detectors (ND) is defined as the distance from 
the Nal(Tl) detector to the neutron detector, while for the 
Nal(Tl) detector it is the distance from the LiF target to the 
Nal(Tl) detector. The neutron detectors are cylindrical Eljen 
detectors with EJ-301 as the scintillator and EJ-510 as the 
reflector. The detector size given is both the diameter and 
length of the cylinder. 


Detector 

Detector 

Scattering 

Recoil 

Flight 


Size (cm) 

Angle (deg) 

Energy (keV) 

Distance (cm) 

Nal(Tl) 

2.5 

0 


50 



0 


91 

ND1 

12.7 

59.1 

29.0 

150 



74.2 

43.0 

135 

ND2 

12.7 

41.3 

15.0 

150 



54.4 

24.9 

122 

ND3 

12.7 

24.9 

5.7 

200 



31.1 

8.8 

164 

ND4 

5.1 

47.9 

19.4 

70 



84.0 

51.8 

52 

ND5 

5.1 

32.2 

9.1 

70 



64.6 

33.3 

41 

ND6 

5.1 

18.2 

2.9 

70 



41.1 

14.3 

33 


Figure [2] Signals from both the Nal(Tl) and the neutron 
detectors were amplified and fanned out. One copy of 
the signals was used to produce a trigger for a CAEN 
V1720E digitizer module (12 bit, 250MS/s), while the 
other copy was digitized. For each trigger, a signal re¬ 
gion of 8 /is with 2 //s before the trigger was digitized 
to ensure that the entire Nal(Tl) scintillation event was 
recorded (scintillation lifetime r=200-300ns). 

The PMT signal from the Nal(Tl) detector was ampli¬ 
fied with a xlO front-end amplifier module developed at 
the Laboratori Nazionali del Gran Sasso (LNGS) while 
the neutron detector signals were sent to a Phillips 779 
xlO amplifier. Amplified signals from both detectors were 
sent to a LeCroy 428F linear fan-out. The signals to be 
used for the trigger first went to low-threshold discrimi¬ 
nators (Phillips 711 and LeCroy 621 AL), whose discrim¬ 
ination levels were set at 1.5 p.e. in order to reach a low 
energy threshold while reducing the random coincidence 
rate. The discriminator outputs for the neutron detectors 
combined by a NIM logical fan-in whose OR output was 
subsequently combined with the Nal(Tl) discriminator 
output in a logical AND (NIM 375L). The coincidence 
window for this AND logic was set to 400 ns to conser¬ 
vatively include the neutron coincidence events with the 
longest time of flight. Subsequent triggers within the ac¬ 
quisition window were discarded. The signal from the 
pulsed proton beam was not used in the trigger, but was 



FIG. 2. Electronics scheme for the measurement. The 
Nal(Tl) signal was fed through a front-end amplifier module 
developed at LNGS. The Nal(Tl) and neutron detectors with 
lower gain were passed through a Phillips 711 discriminator 
with a threshold of 10 mV, while other neutron detectors and 
the proton pulse selector were passed through the LeCroy 
621 AL discriminator with a higher threshold equivalent to 
1-2 photoelectrons in the neutron detectors. 


recorded for off-line analysis. Due to the degradation of 
the proton pulse selector signal in long transmission lines, 
the signal was fed through a discriminator before being 
digitized. 

A basic online analysis was performed to show the scin¬ 
tillation waveforms, the coincidence event rate, and the 
time-of-flight spectra for all neutron detectors. The wave¬ 
form data from all channels were saved to disk in a binary 
format in real time, to be used for offline analysis, as dis¬ 
cussed in Sec. uni 


E. Measurement Summary 

Data were collected in the first position configuration 
for 26 hours and in the second position configuration for 
20 hours, giving approximately 1,000-4,000 coincidence 
events per energy. 

Calibrations of the light yield of the Nal(Tl) detector 
were taken in-run by observing the 57.6 keV gamma ray 
that comes from the first excited state of 127 I, which can 
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be induced through inelastic scattering of the neutrons. 
The light yield of this detector was initially measured to 
be 18.2 p.e./keV ee , but decreased over the course of the 
measurement to 13.7p.e./keV ee due to some degradation 
of the crystal from moisture exposure, and possibly also 
a degradation of the optical coupling. The in-run cal¬ 
ibration compensated for the loss of light yield in the 
calculation of an energy spectrum for the nuclear recoil 
events. 

Separate calibration runs with 133 Ba and 241 Am 
sources also observed this light-yield degradation. How¬ 
ever, a few-percent systematic difference in the light yield 
between the two measurements was observed; the source 
calibrations showed a lower light yield than the real-time 
calibration with 127 I. One potential effect is the skew¬ 
ing of the peak due to the energy of the iodine recoil 
itself, but this effect was estimated to be at or below 1%. 
Another potential reason for this difference stems from 
the position distribution of the scintillation events; the 
gammas from the first excited state of 127 I are evenly 
distributed throughout the crystal, whereas the gammas 
from the external sources will interact within a few mm 
of the crystal edge. This effect can cause a systematic de¬ 
crease in the light yield, as the light yield may be position 
dependent. The in-beam measurement of the 57.6 keV 
peak was used to calibrate our detector performance in 
our analysis. 

After the data were collected for the measurement of 
the quenching factor, a separate run was conducted to 
measure the trigger efficiency of the Nal(Tl) detector at 
low energies. The setup is shown in Figure [3j A Bi¬ 
cron 76-mm Nal(Tl) detector was set up at a reason¬ 
able distance away from the 25-mm Nal(Tl) detector. 
A 22 Na source was placed directly between the two de¬ 
tectors, which produced back-to-back 511-keV gamma 
rays. 133 Ba and 241 Am sources were also put near the 
25-mm detector opposite the 76-mm detector to increase 
the event rate. The 76-mm detector was used as a trigger 
for the measurement with a high threshold. The 25-mm 
detector signal from the discriminator as well as directly 
from the amplifier were fed into the digitizer to measure 
the efficiency of the trigger for low energy recoils. 


III. DATA ANALYSIS 

As discussed in Section [TTJ data from 12 neutron¬ 
scattering angles, corresponding to 12 different nu¬ 
clear recoil energies, were collected in this measure¬ 
ment. The nuclear recoil signals in the Nal(Tl) crys¬ 
tal were selected using time-of-flight (TOF) and pulse- 
shape-discrimination (PSD) cuts. Their energy spectra 
were then compared to the predicted recoil-energy spec¬ 
tra, produced by Monte Carlo simulations, to evaluate 
the energy-dependent nuclear recoil quenching factors. 
We report the analysis of sodium recoil quenching effects 
in an energy window of ^3keV nr to ^52keV nr , which 
covers the DAM A/LIBRA region of interest. Iodine re- 


133 Ba 

O 

241 Am 

O 


25-mm Nal(Tl) Detector 
(used in measurement) 



FIG. 3. Setup for the trigger efficiency measurement of the 
25-mm Nal(Tl) detector. A Na-22 source was placed between 
the Nal(Tl) detector and a separate 76-mm Bicron Nal(Tl) 
detector so that the back-to-back 511keV gamma rays could 
be detected by both detectors. The signal was fed both di¬ 
rectly into the digitizer and also through the discriminator, 
while the Bicron detector was used as a trigger to reduce the 
acquisition of data with no pulses. Additional sources were 
placed near the 25-mm detector to increase the event rate. 


coils were not observed in the measurements due to their 
low recoil energies and larger quenching effects; limits 
were set on the iodine quenching factors. 


A. Data Processing 

The raw data acquired in the measurements contain 
the waveforms of the Nal(Tl) detector, the waveforms of 6 
liquid scintillator neutron detectors, and a periodic pulse- 
selector signal from the proton accelerator, which relates 
to the proton-on-target (POT) time with a constant off¬ 
set. The waveform baselines were first subtracted using 
a drifting-baseline-finding algorithm, which was tuned 
to suppress low-frequency electronic noise while preserv¬ 
ing high-frequency scintillation pulse signals. Individ¬ 
ual pulses with an amplitude >0.2 photoelectrons were 
tagged and further processed for the analysis. For each 
pulse that contributed to a coincidence trigger, all fol¬ 
lowing pulses within 4 /is were clustered together as one 
scintillation event for energy evaluation. The 4 /rs time 
window was chosen to contain the full scintillation signals 
of the highest-energy events considered in this analysis. 

For every coincidence event, we first tried to iden¬ 
tify which neutron detector contributed to the trigger 
by comparing the pulse arrival times with the coinci¬ 
dence trigger time. If the signals from the responsible 
neutron detector and the Nal(Tl) detector satisfied cer¬ 
tain event selection criteria, the Nal(Tl) signal is kept 
for the quenching factor analysis and the neutron detec¬ 
tor position is used to determine the neutron scatter¬ 
ing angle. The most important event selection criteria 
in the analysis is the cut on the time of flight (TOF), 
or the time difference between the pulse arrival times 
in different detector channels. For the neutron energy 
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T0F1 (ns) Energy (keVee) 

FIG. 4. Left: The double TOF spectrum for coincidence events containing ~29keV nr Na recoils. Gammas (first vertical 
band) and neutrons (second vertical band) from the LiF can be separated by TOF1, defined as the TOF from LiF to Nal(Tl); 
neutron scattering off Nal(Tl) can be further selected (blue box) using TOF2, defined as the total TOF from the LiF to the 
liquid scintillator neutron detectors. Right: The energy distribution of the neutron induced events from the second band in 
the figure on the left. The blue box contain the neutron scattering events. Neutron-induced nuclear recoils are marked in the 
oval and the 127 1 excitation events are shown in the vertical band. 


used in this measurement (^690keV, corresponding to 
a speed of ~ 1 cm/ns), the time required for the neu¬ 
trons to travel from the Nal(Tl) detector to the neutron 
detectors was at the level of ^50 -200 ns, depending on 
the detector positions. This well-defined time correla¬ 
tion made the neutron events distinct from the nearly 
instantaneous gamma coincidence background and ran¬ 
dom coincidence backgrounds. Furthermore, the pulsed 
neutron beam also made it possible to calculate and cut 
on the TOF between the LiF target and the Nal(Tl) de¬ 
tector, which further suppressed the gamma ray back¬ 
ground generated by the proton beam. 

In this analysis, TOF1 was defined as the difference 
between the arrival times of the Nal(Tl) signal and the 
proton pulse signal, and TOF2, or the total TOF, was 
defined as the difference between the arrival times of the 
neutron detector signal and the proton pulse signal. In 
a real neutron-induced coincidence event, TOF1 repre¬ 
sents the time (with a constant offset) required for the 
neutron to travel from the LiF target to the Nal(Tl) de¬ 
tector, and TOF2 corresponds to the time (with a similar 
constant offset) required for the neutron to travel from 
the LiF target, to scatter off the Nal(Tl) detector, and 
to be recorded by a liquid scintillator neutron detector. 

Figure [4] (left) shows the double TOF spectrum 
(TOF2 vs. TOF1) for the coincidence events between 
the Nal(Tl) detector and a neutron detector contain¬ 
ing ^29keV nr Na recoils. The first and second verti¬ 
cal bands, respectively, correspond to the gammas and 
neutrons that were produced by the proton beam and 
recorded by the Nal(Tl) detector. Because the gamma- 
ray flight time from the target to the Nal(Tl) detector 
is only a few nanoseconds, the time separation between 
the two vertical bands provided an estimate of the neu¬ 
tron TOF, which was confirmed by direct calculations. 
The experiment was designed in a way that the gamma 


and neutron TOF bands were sufficiently separated, so 
conservative TOF cuts could be used to efficiently reject 
the gamma background with little impact on the neu¬ 
tron events. Events in the horizontal band with TOF2 
slightly below -0.15/is were identified to be the beam- 
induced gamma rays that directly hit the neutron detec¬ 
tor in coincidence with a random Nal(Tl) scintillation 
event. The diagonal band with TOF2^TOFl was at¬ 
tributed to simultaneous scintillations in the Nal(Tl) de¬ 
tector and in the neutron detector from environmental 
radioactivity such as high-energy gammas or comic-ray 
showers. Combining the analysis using both TOF1 and 
TOF2, the blue box in Figure [4] (left) was identified to 
contain the desired neutron coincidence events where a 
neutron scattered off the Nal(Tl) detector and then got 
recorded in the neutron detector. 

In addition to nuclear recoils, neutron interactions with 
Nal(Tl) also produced nuclear excitations via inelastic 
scattering. As shown in Figure [4] (right), the 57.6keV 
127 1 excitation gamma rays were observed (vertical band 
in the plot) in all neutron induced Nal(Tl) scintillation 
energy spectra. As discussed in Section [TlJ these gam¬ 
mas were used to provide an in-run energy calibration 
for this measurement; they also provided a way to moni¬ 
tor and correct the degradation of the Nal(Tl) light yield 
observed between the runs. It was estimated that this 
calibration introduced ~ 1 % uncertainty in the energy 
scale due to the iodine recoils accompanying the 57.6 keV 
gamma rays, and the time-dependent light yield correc¬ 
tion further introduced a 1.5% uncertainty. The quench¬ 
ing factors to be reported in this paper are all normalized 
to the scintillation efficiency of Nal(Tl) under 57.6 keV 
gamma excitations. Due to a possible non-linearity of 
the Nal(Tl) scintillation output |27j 28], the evaluated 
nuclear-recoil quenching-factor values may depend on the 
gamma calibration point, and the results from different 
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measurements need to be appropriately scaled for direct 
comparison. 

With the neutron events selected using TOF cuts and 
the energy scale calibrated with 127 I excitation gammas, 
the nuclear recoil energy spectra were extracted for all 12 
neutron scattering angles. The single-scattering Na re¬ 
coil peaks could be resolved clearly in the relatively high- 
energy regions (>10keV nr ), as illustrated in Figure [5] 
(left), which shows the energy spectrum of ~ 29keV nr Na 
recoil events selected from Figure [4] However, the energy 
spectra of relatively low-energy Na recoils (<10keV nr ) 
were contaminated by noise from random coincidence 
backgrounds, as shown in Figure [5] (right). Fortunately, 
liquid scintillators have a scintillation pulse shape that 
depends on the particle type [27], a fact that we ex¬ 
ploited to reject gamma backgrounds. As shown in the 
insert figure of Figure [5] (right), a conservative cut on 
F50, defined as the fraction of the pulse integral in the 
first 50 ns, was found to be sufficient in rejecting the low 
energy background events. A similar F50 cut was applied 
to the Nal(Tl) detector signals to reject fast Cerenkov 
light in the Nal(Tl) crystal and fast noise signals in the 
electronics; typical Nal(Tl) scintillations have a lifetime 
of ^200-250 ns and the bias from this conservative F50 
cut was negligible. We comment that the PSD cuts were 
only necessary for the analysis at the lowest energies 
(<1.5keV ee ), below the DAMA/LIBRA energy thresh¬ 
old of 2 keV ee . 

With the double-TOF cuts and the PSD cuts at low 
energies, the spectral peaks from single neutron-Na scat¬ 
tering were identified at all neutron scattering angles ex¬ 
cept for the 2.9keV nr recoils. We comment that all cuts 
used in this analysis were carefully chosen to be overly 
conservative so the analysis of neutron-induced events 
should not be biased. Thanks to the application of mul¬ 
tiple event-selection criteria including double-TOF cuts 
and double-PSD cuts, the conservative cuts, when com¬ 
bined, were demonstrated to be sufficiently efficient in 
suppressing background events that were not due to neu¬ 
tron scattering coincidences. 


B. Quenching Factor Evaluation 

To evaluate the Na recoil quenching factors, the 
recorded nuclear recoil spectra need to be compared to 
the expected recoil spectra. In principal this is straight¬ 
forward to calculate with the given kinematics. However, 
a number of effects could introduce uncertainties, includ¬ 
ing the proton energy loss in the LiF target, the energy 
dependence of the 7 Li(p,n) 7 Be reaction, the finite size 
of the Nal(Tl) detector, and the finite size and angular 
positions of the liquid scintillator neutron detectors. In 
this analysis, we employed Monte Carlo techniques to cal¬ 
culate the energy spectra of sodium and iodine nuclear 
recoils, taking into account these effects. 

The Geant4 package [29] (version 4.9.6.p3) was used 
to simulate the neutron interactions in the experimental 


setup. The beam-produced neutrons were generated us¬ 
ing the 7 Li(p,n) 7 Be cross-section data measured by Burke 
et al. [26] . with the proton energy loss in the LiF target 
appropriately simulated. Thanks to the small size of the 
Nal(Tl) crystal used in this measurement (a cube with 
25-mm sides), single-scattering nuclear recoils dominated 
the simulated energy spectrum, as illustrated in the insert 
of Figure [6] while multiple-scattering events contributed 
a small, featureless, near-flat background. 

The Na-recoil quenching factors were extracted by fit¬ 
ting the observed nuclear recoil spectra to the Monte 
Carlo-simulated energy spectra convolved with a Gaus¬ 
sian spread. The maximum-likelihood algorithm was 
used to appropriately handle the fit in the low statis¬ 
tic data bins, and all fits returned reasonable y 2 values. 
An example of the spectral fit is shown in Figure [6] along 
with the simulated energy spectrum at the corresponding 
neutron detector configuration (the insert of Figure [6]). 
The systematic uncertainty in the fit process was evalu¬ 
ated by varying the fit range and comparing the fit re¬ 
sults. 

Based on the trigger-efficiency measurement described 
in Section [II] we found that the Na-recoil events at low 
energies (< 10keV nr ) were subject to loss of triggers due 
to the small pulse heights falling below the discriminator 
threshold. As shown in Figure [7] the trigger efficiency 
drops dramatically around 20 - 25 digitizer counts, cor¬ 
responding to ^1.5 photoelectron amplitude. It was es¬ 
timated that approximately 50% of the Nal(Tl) scintil¬ 
lation events at 0.65 keV ee were lost due to the threshold 
effect and the loss was more significant at lower energies. 
This threshold effect slightly biased the low-energy tail 
of the 9.1 keV nr Na recoil data, as illustrated in Figure [7] 
but impacted the 8.8keV nr , 5.7keV nr and 2.9keV nr Na 
recoil data to a larger extent. 

Therefore, we adopted a different analysis strategy 
at these energies. First, the measured trigger-efficiency 
curve was used to restore the loss of triggers for events 
with > 10% trigger efficiency. If the peak of the pulse- 
height spectrum could be resolved after this correction 
(5.7keV nr , 8.8keV nr and 9.1keV nr ), the peak position of 
the pulse-height spectrum was used to estimate the peak 
position of the pulse-integral (energy) spectrum. The 
pulse height/integral correlation was obtained by inves¬ 
tigating the pulse-energy distribution at different pulse- 
height values. This method was confirmed to yield the 
correct peak energy in tests with higher-energy recoils 
that did not suffer trigger loss. The quenching factors 
were then calculated by comparing the peak positions of 
the observed energy with that of the predicted energy. 
For the Na recoils of the lowest energy (2.9keV nr ), the 
pulse-height spectrum could not be effectively restored, 
so an upper limit of the quenched energy (and the corre¬ 
sponding quenching factor) was extracted. 

The Na-recoil quenching results are summarized in Ta¬ 
ble |H[ Due to an asymmetry in the energy spectra, the 
fitted quenching factor values differ slightly from the esti- 
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FIG. 5. Examples of observed Na recoil energy spectra after the conservative TOF cuts. Left: Energy spectrum of ~29keV nr 
Na recoils (selected from Figure [4|. Only cuts on TOF1 and TOF2 were applied. Right: The energy spectrum of ~5.7keV nr 
Na recoils with (solid blue) and without (dotted red) PSD cuts on the liquid scintillator detector signals. The insert figure 
shows the PSD parameter F50 vs. the pulse integral (in the unit of number of photoelectrons, or NPE) in the neutron detector. 
Note that PSD cuts were only necessary for data below the DAMA/LIBRA energy threshold of 2keV ee - 



FIG. 6. Spectral fit of ~15.0keV n r Na recoil events with 
the Monte Carlo-simulated spectrum with a Gaussian spread. 
The measured recoil spectrum is shown in black and the fit 
function is shown in red. Insert: The Monte Carlo-simulated 
nuclear recoil energy spectra for the corresponding neutron 
scattering angle. The blue spectrum shows the energy of 
single-scattering sodium recoil events; the red also includes 
iodine recoils and multiple scattering events in Nal(Tl). The 
background rate of iodine recoils and multiple scatterings is 
~ 20 times lower than that of single-scattering sodium recoils 
in the relevant energy region. The peak around 3 keV nr in the 
insert figure corresponds to single-scattering iodine recoils, 
which fell below the trigger threshold in the measurement. 
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FIG. 7. The pulse-height spectra of 5.7keV nr Na recoil events 
(top) and 9.1 keV nr Na recoil events (bottom). The measured 
pulse-height spectra are shown in blue, and the independently 
measured trigger efficiency is shown in black dotted line (ef¬ 
ficiency scale is shown with the axis on the right). 


mates based on the peak positions of the spectra, and this 
uncertainty was included in the peak-comparison analy¬ 
sis at low energies. In addition to the uncertainties from 
the spectral fits and peak-comparison analysis, the re¬ 
sults also include the 1.5% uncertainty from the gamma- 
ray calibration, and 3-12% uncertainty from the detector 
position measurements, which varied with neutron scat¬ 
tering angles and distances between detectors. 

Although this measurement was designed to study the 


quenching effect of low-energy sodium recoils, data ac¬ 
quired with large scattering angles could also contain io¬ 
dine recoils of up to 10keV nr based on kinematic calcu¬ 
lations and simulations. However, we did not observe 
significant evidence for iodine recoils with the expected 
rate above 0.65keV ee , in which region we have over 50% 
trigger efficiencies. Therefore, we have set an upper 
limit of 0.065 (>3cr) for the iodine quenching factor at 
10keV nr . We note that DAMA uses iodine quenching 
value of 0.09 PH. and Collar measured a much lower 
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TABLE II. Summary of the Na quenching factors measured 
in this work. The angles are calculated using the central po¬ 
sitions of the detectors; the energies reported are the peak 
values and widths. Quenching factors were evaluated by spec¬ 
tral fits between observation and simulation above 10keV n r 
and by comparing the peak energy positions at lower energies 
after correcting for the trigger-efficiency loss, as described in 
the text. 


Scattering 
angle (°) 

Sim. Na recoil 
energy (keV nr ) 

Observed recoil 
energy (keV ee ) 

Quenching 

factor 

18.2 

2.9 ± 0.7 

<0.65 

< 0.22 

24.9 

5.7 ± 0.7 

0.76 ± 0.4 

0.133 ± 0.018 

31.1 

8.8 ± 1.2 

1.13 ± 0.5 

0.129 ± 0.014 

32.2 

9.1 ± 1.2 

1.46 ± 0.5 

0.162 ± 0.012 

41.1 

14.3 ± 2.4 

2.21 ± 0.9 

0.159 ± 0.019 

41.3 

15.0 ± 1.4 

2.36 ± 0.8 

0.160 ± 0.010 

47.9 

19.4 ± 1.6 

3.21 ± 1.0 

0.168 ± 0.009 

54.4 

24.9 ± 2.4 

4.10 ± 1.5 

0.171 ± 0.010 

59.1 

29.0 ± 1.9 

5.36 ± 1.9 

0.188 ± 0.008 

64.6 

33.3 ± 2.8 

6.19 ± 2.1 

0.191 ± 0.011 

74.2 

43.0 ± 2.2 

8.53 ± 2.7 

0.204 ± 0.008 

84.0 

51.8 ± 2.6 

10.59 ± 4.5 

0.207 ± 0.010 
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FIG. 8. The Na-recoil quenching factors measured in this 
work in comparison with previous results mm- With much 
improved accuracy, this measurement approximately agrees 
with previous measurements at high energies and may recon¬ 
cile the low-energy discrepancy. We comment that these mea¬ 
surements were calibrated to different gamma ray energies, 
and this direct comparison of quenching factors may include 
an uncertainty from the possible non-linearity of Nal(Tl) scin¬ 
tillation nail]. 
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IV. RESULTS AND DISCUSSIONS 

By using low-energy, pulsed neutrons and combining 
TOF and PSD methods, we obtained the most accurate 
measurement of Na-recoil quenching factors in Nal(Tl) to 
date. This result differs significantly from the measure¬ 
ments of DAMA, as illustrated in Figure [8] At high en¬ 
ergies, our measurement approximately agrees with pre¬ 
vious reports, but the quenching factor values are found 
to drop significantly as the recoil energy decreases, simi¬ 
lar to what was observed in the Collar measurement [18] . 
At low energies, our measurement falls between, and is 
statistically consistent with, the recent measurements by 
Collar and Chagani m ■ This energy-dependent quench¬ 
ing effect may be explained by the increasing ionization 
density along the Na recoil tracks at low energies, which 
has a large nuclear component but small electronic com¬ 
ponent; it is usually believed that the former leads to 
heat generation and the latter causes scintillation na. 

The Na-recoil quenching factor decreases from ~0.19 
at 6keV ee to ^0.15 at 2keV ee , which means that the 
DAM A/LIBRA modulation signal (most notable at 2 
- 6keV ee ), if attributed to Na recoils, occurs in the 
energy window of 13 keV nr - 32 keV ni 0 instead of the 


1 Here we have assumed that our quenching factors calibrated 


7 keV nr - 20 keV nr window, as was previously thought. 
Due to the exponentially decreasing characteristic of a 
WIMP scattering spectrum, the WIMP explanation of 
the DAMA/LIBRA signal would require the WIMPs 
to have larger masses and/or larger interaction cross 
sections with nucleons than those assumed in previous 
DAMA/LIBRA interpretations. 

Using a standard WIMP-halo model and considering 
only spin-independent WIMP-nucleon interactions EU 
31 , we evaluated how this new quenching measure¬ 
ment would impact the dark-matter interpretation of 
the DAMA/LIBRA signal. Figure [9] shows the WIMP 
fits of the DAMA/LIBRA modulation data [1! with the 
new quenching factor values, compared with those us¬ 
ing old quenching values. Since this work focuses on the 
Na quenching effects below 52 keV nr (10.6 keV ee ), we re¬ 
stricted the fits to below 10keV ee energy, above which 
the DAMA/LIBRA signal is consistent with zero modu¬ 
lation. As expected, the WIMP analysis yielded a global 
best fit with a heavy WIMP (^70GeV/c 2 mass) and a 
local best fit with a light WIMP 10 GeV/c 2 mass), as 
shown in Figure [l0| 

With the new quenching factors, the light-WIMP fit 
(dominated by Na recoils) does not represent the data 
well, producing a Xmin °f 36 with 14 degrees of freedom 
(P<0.01) in comparison with a X^in °f 19 (P~0.2) for 


to the 57.6 keV gammas can be directly applied to the 
DAMA/LIBRA data because DAMA/LIBRA uses gamma rays 
of a similar energy (59.5 keV) for regular calibrations. 










































10 


Oh 





Energy (keV ) 


FIG. 9. Spectral fits of the DAMA/LIBRA modula¬ 
tion amplitudes to the standard WIMP model (spin- 
independent only) with the new Na quenching factors in 
comparison with the fits with the old values. The heavy- 
WIMP fit does not change, but the low-mass-WIMP pic¬ 
ture no longer fits the data well. 
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FIG. 10. The DAMA/LIBRA la, 3a, 5 a significance 
contours in the WIMP parameter space. The newly- 
obtained results using the quenching factors measured 
in this work are shown in color-filled regions, while the 
old results with the DAMA/LIBRA quenching factors 
are shown in colored lines. The heavy-WIMP contours 
do not have significant change, but the 1 a and 3 a con¬ 
tours in the low-mass WIMP regions disappeared com¬ 
pletely, disfavoring a light WIMP. The dashed line is the 
dark matter exclusion curve calculated using the overall 
DAMA/LIBRA observed event rate J22]- In the standard 
WIMP picture, WIMP parameters above this line would 
produce a nuclear recoil spectrum above the observed one 
in DAMA/LIBRA at > 1 energy bin. 


the standard DAMA/LIBRA fit with the old quench¬ 
ing value. Therefore, the low-mass WIMP region is 
strongly disfavored, as illustrated by the diminishing 


10 


X 2 significance contours around 10 GeV/c 2 in Figure 
Moreover, by shifting the light-WIMP contours to larger 
WIMP-mass and cross-section values, the tension be¬ 
tween DAMA/LIBRA and other experiments increases 
in the standard WIMP picture. 


The high-mass-WIMP fit does not change significantly 


because the recoil signal is dominated by WIMP-iodine 
scatterings and the effect due to sodium quenching is 
negligible. Nonetheless, the best fit values in the high- 
mass WIMP region lead to a WIMP-interaction rate in 
Nal(Tl) higher than that observed in DAMA/LIBRA 
around 2 keV ee [3T . The dark matter exclusion curve m 
in Figure [To] (dashed line) was calculated using the ob¬ 
served event rate in DAMA/LIBRA between 2keV ee 
and 6 keV ee , and in the standard WIMP picture, any 
WIMP parameters above this line would produce a nu¬ 
clear recoil event rate higher than what was observed 
in DAMA/LIBRA at one or more energy bins. More 
importantly, the iodine-dominated heavy-WIMP region 
has been excluded independently of WIMP models by 
the KIMS experiment using CsI(Tl) crystals [33j. 

However, the standard WIMP models are known to 
have large uncertainties, and as discussed earlier, alter¬ 
native dark matter theories may still be able to reconcile 
the experimental results ia no], a model-independent 
test of DAMA/LIBRA, therefore, is best made with a 
Nal(Tl) experiment with lower background than that of 
DAMA/LIBRA. 

The lattice orientation of the Nal(Tl) crystal used in 
this experiment was not measured, so this measurement 
does not provide a sensitive test of the possible ion¬ 
channeling effect [34]. But, with the large number of 
Na-recoil angles measured and the fact that all quench¬ 
ing factors line up on a curve well below unit quench¬ 
ing (no quenching), we do not observe any evidence for 
the channeling effect. Similarly, although the setup for 
this experiment was optimized to measure sodium recoils, 
iodine recoils should have been observed if the quench¬ 
ing factor were at the value measured by DAMA-LIBRA 
(0.09). Based on the absence of the iodine recoil peaks, 
we set an upper limit (>3cr) of 0.065 on the iodine recoil 
quenching factor at 10keV nr . 


V. CONCLUSION 

We carried out an accurate measurement of the rela¬ 
tive Nal(Tl) scintillation efficiency for Na recoils induced 
by a pulsed neutron beam, covering an energy window 
from 3keV nr to 52keV nr (or 0.65 - 10.6keV ee , covering 
the whole DAMA/LIBRA modulation signal region of 
2 - 6 keV ee ). By using double-TOF cuts and double- 
PSD cuts, we suppressed the coincidence background in 
the measurement with a high efficiency and obtained the 
most accurate results to date. 

The Na recoil quenching factors are found to de¬ 
crease significantly at low energies, which caused the 
DAMA/LIBRA modulation signal to be less compatible 
with a light-WIMP explanation in the standard WIMP 
picture. Although alternative models may still be able 
to reconcile the DAMA/LIBRA signal with other experi¬ 
mental results, a model-independent test using ultra-high 
purity Nal(Tl) crystal detector is necessary to confirm or 
refute the DAMA/LIBRA dark matter claim. 
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